Effects of static and temporally fluctuating tensions on semiflexible polymer looping 
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Biopolymer looping is a dynamic process that occurs ubiquitously in cells for gene regulation, 
protein folding, etc. In cellular environments, biopolymers are often subject to tensions which 
are either static, or temporally fluctuating far away from equilibrium. We study the dynamics of 
semiflexible polymer looping in the presence of such tensions by using Brownian dynamics simulation 
combined with an analytical theory. We show a minute tension dramatically changes the looping 
time, especially for long chains. Considering a dichotomically flipping noise as a simple example of 
the nonequilibrium tension, we find the phenomenon of resonant activation, where the looping time 
can be the minimum at an optimal flipping time. We discuss our results in connection with recent 
experiments. 



I. INTRODUCTION 

The kinetics of polymer loop formation has been stud- 
ied for several decades and recently has attracted re- 
newed attention due to the particular importance in bi- 
ology. The DNA loop formation is a basic process that 
underliesgenetic expression, replication, and recombina- 
tion [l|, For example, in E. coli the lac repressor 
(Lacl)-mediated loop is crucial for the repressive regu- 
lation of lac genes. The hairpin loop formation is the 
elementary step in protein folding 0] and structure for- 
mation in RNA folding 

A cell is crowded with a multitude of subcellular struc- 
tures including globular proteins and RNAs with 
which DNA is constantly interacting. A DNA fragment 
about to loop is often subject to temporally fluctuating 
forces due to its dynamic environment including the other 
part of the chain. Recently, the power spectrum of the 
fluctuating force exerted on cytoskcleton was measured to 
be an order of magnitude larger than that expected from 
thermal equilibrium condition This indicates that 
the cell interior is an active and nonequilibrium medium. 

The advance of single molecule experiment techniques 
provides detailed information on the DNA loop forma- 
tion. Finzi and Gelles Q observed Lacl-mediated DNA 
loop formation and dissociation by monitoring nano-scale 
Brownian motion of the micron-sized particle attached to 
one end of the DNA. Lia et al. Q showed that in gal re- 
pressor and DNA-bending protein HU mediated looping, 
mechanical constraints such as tension and torsion play a 
pivotal role. Gemmen et al. [§| studied effects of tension 
in the presence of two-site restriction enzymes which can 
cut the DNA upon binding on two sites simultaneously 
They found that the cleavage activity decreases approx- 
imately 10-fold as the tension increases from 30 fN to 
700 fN. They also found that the optimum loop size de- 
creases with the tension, which is qualitatively in agree- 
ment with theoretical predictions jiol ]. More recently, 
Chen et al. [llj studied effects of tension in femtonew- 
ton range on the kinetics of Lacl-mediated DNA looping. 
They found that small tension of 100 fN scale on the sub- 
strate DNA can not only increases the looping time [ll[ 
but also found that the looping time is greatly reduced in 



the presence of a fluctuating tension [12|. These results 
suggest the ubiquitous roles of the static and temporally 
fluctuating tensions in regulation of the DNA loop for- 
mation. Yet, there appears to be no unifying conceptual 
or theoretical framework that explains a variety of exper- 
iments including these. 

Theoretically, on the other hand, Yan et al. [l3[ devel- 
oped a transfer matrix method to calculate semiflexible 
polymer end-to-end distance distribution function and 
loop formation probability (or J- factor). They studied 
various effects of nonlinear elasticity arising from DNA 
bending protein-induced kinks [H or thermal-fluctuation- 
induced bubbles on the J-factor. Their study provides a 
valuable insight to understand DNA bending on short 
length sca le 1 131 , which has attracted much attention re- 
cently ll^lTa . They also studied effects of tension on the 
J-factor [13( | . which is related to the free energy barrier 
for loop formation thus to the loop formation rate. 
Similar results are obtained by using an elastic theory 
of a semiflexible polymer [l(|. However, since the loop 
formation rate is not proportional to the J-factor alone 
but depends on the free energy given an arbitrary chain 
end-to-end distance, it is hard to quantitatively compare 
these theories to the experiment [11|. 

Independently, Blumberg et al. [l7| studied effects of 
static tension on protein-mediated DNA loop formation 
by modeling the DNA conformation to be either one of 
two states, looped and unloopcd states. In their appeal- 
ing calculations of free energy change associated with the 
transition under the tension, they considered not only 
the stretching free energy of DNA but also DNA align- 
ment constraint imposed by protein binding. They found 
that for the loop size larger than 100 base pair distance 
(w 34 nm), a tension of 0.5 pN can increase the loop- 
ing time by more than two order of magnitude. There is 
room for improvement in their approach, however, on the 
evaluation of the free energy that can be valid for short 
end-to-end distance of the chain as well as a description 
of detailed kinetic process using the mean first-passage 
time approach. 

In this paper, in an effort to understand the basic 
physical mechanism of the biopolymer looping in a coher- 
ent manner, we perform Brownian dynamics simulation 
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of semiflexible polymers treated as extensible wormlikc 
chain, combined with one-dimensional theory of barrier 
crossing over the free energy of loop formation. For ana- 
lytical understanding, we use, as an example, the mean- 
field wormlike chain model jig , which is shown to be 
a good approximation for the free energy for the chain 
lengths we consider here. With static tensions, we find 
that the looping time, defined as the mean first-passage 
time to cross the free energy barrier, steeply increases 
with the applied tension /, in an agreement with our 
simulation results but distinct from the previous theoret- 
ical result [13 • For the case of time-dependent tension, 
we consider dichotomically fluctuating tension, where the 
looping times are found to be reduced, consistent with 
the experiment [Hj]. Most importantly, we find so-called 
the resonant activation, where the looping time is the 
minimum at an optimal flipping time of the dichotomic 
force. In this exploratory study, we neglect the alignment 
constraint on the loop formation, which is minor effect 
for the chain lengths we consider here [lTj ■ 

In the following section, we describe our polymer model 
and simulation method, whose results are discussed in 
Sec. IIIII We summarize our results in Sec. IIV1 



II. POLYMER MODEL AND SIMULATION 
METHOD 

We consider the semiflexible polymer looping in the 
presence of static and fluctuating tension. The polymers 
are modeled as semiflexible chains of N beads of diameter 
d, with the interaction potential U. Here U = U s + Ub, 
where U s and Ub are the stretching and bending energy 

k N ~ 1 

U s = ^ £(|* + i -Fil-Zo) 2 , (1) 



u b = ^£el (2) 

i=2 

where k is the stretching stiffness, Iq is the natural bond 
length, k is the bending stiffness, r*j is the position of 
the ith bead, and Oi is the angle of the ith bond. The 
dynamics of N beads (i = 1,2, ...,7V) are described by 
the overdamped Langevin equation 

7f$(t) = -ViU + £(t) + f{t)(S hN - 6i,i), (3) 

where 7 is the friction coefficient of the bead, and £j is 
the Gaussian random noise with mean (£»(t)) = and 
variance £7,9(0)} = 2jk B T5ij6p. q 6(t). Here, (...) 

denotes ensemble average, p and q are the Cartesian co- 
ordinate indices, k B is the Boltzmann constant, and T 
is the absolute temperature. The additional force f(t) 
is the tension which is applied only to the two end seg- 
ments (i = 1 and N) along the chain end-to-end direc- 




FIG. 1: A schematic picture of the model semiflexible polymer 
under a negative tension f(t). 



tion, f(t) = f(t)r^- (seeFig.p. The bending stiff- 

\r N - n| 

ness of the short semiflexible chains we consider allows us 
to neglect the excluded volume effect. Furthermore, we 
neglect the hydrodynamic interactions, which is found to 
be small in the short chains we consider here fl9| . 

We use the parameters Iq, k B T, and ^^/(ksT) to fix 
the length, energy, and time scales, respectively. The 
dimensionless parameters in our simulation are Iq = 1, 
k B T = 1, k = 100, and k = 5. In our model, we consider 
lo as 10 nm, so with k = 5 the persistence length l p = 
nlo/(k B T) of the chain is 50 nm. The diameter of the 
bead is set to be 5 nm to fit the hydrodynamic friction 
coefficient of a cylinder with 2 nm diameter and 10 nm 
length, which is about 4.61 x 10 -11 kg s _1 with the water 
viscosity 77 = 0.001 kg m _1 s _1 at room temperature. 
The time unit is then 1.14 /is. The equations of motion 
are integrated by using a second-order stochastic Runge- 
Kutta algorithm [2(j with the time-step At = 2 x 10~ 3 . 
In our simulation, we consider that the looping occurs 
whenever the chain end-to-end distance r = |nv — r±\ is 
shorter than a cutoff distance l c . The average were taken 
over at least 2,000 and 5,000 independent runs for static 
and fluctuating tension cases, respectively. 

III. RESULTS AND DISCUSSION 

A. Effects of static tension on the looping 

We first study the effects of static tension on semiflcx- 
ible polymer looping. Here we consider the bead number 
N = 12, 18, and 24 which respectively correspond to 
the chain lengths L = 110, 170, and 230 nm, similar to 
the loop size considered in recent experiments [ill, [13 ■ 
Figure [5] shows the looping time T as a function of ten- 
sion / for the chain of JV = 12 and l c = 2Zo- We apply 
the tension / from -81 to 162 fN, comparable in magni- 
tude to the typical entropic force on the double-stranded 
DNA, k B T/l p w 80 fN. The / is much smaller than the 
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FIG. 2: The looping time T(/) as a function of tension / for 
the chain of N = 12. 



piconewton scale forces those are typically involved in ac- 
tive molecules 



2l| e.g., molecular motors in a cell. As 



/ changes over the scale as small as 100 fN, the looping 
time T dramatically changes; for example, 7~ increases 
about 5 times as the tension / increases 120 fN. 

To understand this behavior, we consider scmiflcxiblc 
polymer looping as a one-dimensional barrier crossing 
process [lj|, which is described by the Langevin equa- 
tion 



(4) 



Here, F(r) = — k^T log P(r) is the free energy of the 
chain given the chain end-to-end distance r, where P(r) 
is the radial distribution function, and is a random 
force due to thermal fluctuation given by a Gaussian and 
white noise that satisfies = and (£(f)£(0)) = 

'yksTSit). The shape of F{r) obtained by the iV-beads 
simulation is shown in Fig. [3] (solid line with circles) for 
the chain of N = 12. In this one-dimensional description, 
the looping time is the mean first-passage time (MFPT) 
for the variable r to reach the cutoff distance l c starting 
from the initial chain end-to-end distance tq. It is given 
by |H 



drcxp{(3F(r))- J dr' exp(- ( 9F(r')), (5) 

where D = 2ksT/j is the relative diffusion 
coefficient of two end beads, j3 = (fe^T) -1 . 
Then the looping time is T — (T{ro)) eq = 
J dr exp(-(3F(r ))T{r )/ J dr exp(-/3F(r )), where 
(...) q represents the average over the initial equilibrium 
distribution. The MFPT in the presence of tension 
/, also can be calculated by using the free energy, 




FIG. 3: The free energy -Fo(r) of -/V = 12 chain obtained from 
the simulation (solid line with circles) and from the mean-field 
wormlike chain model [ll| (solid line). 



F(r) = Fq(t) — fr, where Fo(r) is the free energy 
without tension. 

For the P(r), we use the mean- field wormlike chain 
model (MF-WLC) [l8| as an example. The radial dis- 
tribution function P(r) of this model is P(r) ~ r 2 [l — 



L> J ' GX IH 4i p (l-(r/L) 2 

sonablc approximation to our simulation result for Fq{t), 
except for r L, as shown in Fig. [3] The large devia- 
tion between two curves in the region near r/L = 1 is 
because our model allows the chain extensibility while 
MF-WLC does not. However, according to the single 
molecule experiment using optical tweezers [23[ , the DNA 
is extensible, manifesting overstretching transition sub- 
ject to a strong tension, thus showing a significant devia- 
tion from the incxtensible WLC. The optimal value k to 
fit the force-extension curve is found to be k « 2700 [23| . 
However, the simulation using k as large as k = 2700 
demands computing time much longer than that using 
k = 100. Unlike the free energy, the normalized looping 
time 7Ar(/)/T7v(0), the looping time in the presence of 
tension / relative to the looping time in its absence, is 
expected to be quite insensitive to the value of k [24| HH . 
For this reason we adopt k = 100, which was also em- 
ployed in a number of studies on the looping [H, [2(| . 

To understand the effect of tension on the looping time, 
the MF-WLC provides a useful analytical model. Fig- 
ureHJshows 7iv(/)/7/v(0) as a function of / for the chain 
lengths iV=12, 18, and 24. Here the symbols are simu- 
lation results and the dashed lines are from MFPT cal- 
culations (Eq. [S]). For the range of / we study here, two 
results are in an excellent agreement. This figure also 
shows that the normalized looping time increases expo- 
nentially with /. The reason is that, for the short chain 
lengths considered here, the free energy of the loop forma- 
tion increases with /(r) oq « fL, so T increases approx- 
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r]. This formula yields a rea- 
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FIG. 4: The normalized looping time 7jv(/)/7jv(0) as a func- 
tion of tension / for N — 12, 18 and 24. The symbols are sim- 
ulation data and the dashed lines are from mean first-passage 
time calculations using mean- field wormlike chain model [lg| . 
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FIG. 5: The relative looping time T(/)/T(/ = 60 fN) as a 
function of tension / with the cutoff distance l c = 5 nm for 
L = 103.7 nm obtained by the MFPT calculation using MF- 
WLC model [ll| (solid line). The circles are experimental 
data including error bar from Ref. pU |. 



imately in exponential with /. In contrast, the theory 
of Ref. [l7| considered that, for low tension (/ < 80 fN), 
7" increases exponentially with f 2 because of Gaussian 
force-extension relation they used. 

The looping time in the presence of a tension / can be 
written in dimensionless form 

where x = r/L. ip is a function of a dimensionless scaling 
variable (3fL. It means for very large L, a minute ten- 
sion / can dramatically change the looping time. Indeed, 
Fig. @] shows that, for longer chains, the normalized loop- 
ing time changes more sensitively with /. If N is 100, 
corresponding to L « 1 fj,m, a change of the force as 
small as 4 fN can affect the looping time. This is a con- 
sequence of the cooperativity of the long polymer chains 
arising from the chain connectivity, which is previously 
addressed in a study of polymer translocation through 
membrane [27| . The sensitivity of looping on tension is 
an emergent behavior that manifests beyond the com- 
plexity of real DNA loop formation (28|, e.g., the details 
of chain conformations outside of the loop, orientation 
constraint, and associated proteins. 

Finally, we compare our result with the recent exper- 
imental data [ll[ for the DNA loop size L = 103.7 nm. 
As shown by the circles in Fig. [5J the looping time rises 
steeply as a function of tension / relative to the one with 
/ = 60 fN, which is the minimum tension used in the ex- 
periment. We also plot the corresponding relative MFPT 
by a solid curve. With the cutoff distance chosen to 
be the bead diameter, l c =5 nm, our result is in good 
agreement with the experiment [111 ]. Because our model 



does not consider many details of DNA loop formation, 
this agreement could be somehow fortuitous but encour- 
aging. A theory including effectively the complexity of 
DNA loop formation and further controlled experiments 
of various loop sizes are needed for quantitative compar- 
ison. 

v(PfL, l j-, l j-U6) 

B. Effects of dichotomically fluctuating tension on 
the looping 

Now suppose that tension temporally fluctuates due 
to the nonequilibrium noise inherent in vivo systems 
which can be generated by a variety of constituents of 
a cell, e.g., protein like RNA polymerase. Also a DNA 
fragment about to loop is influenced by the other part 
of the chain whose conformation is constantly fluctu- 
ating. As a simple example of nonequilibrium fluctua- 
tions, we consider a dichotomic noise, with which the 
tension f(t) flips between two level of the forces +fd 
and —fd with a mean flipping time r. The f(t) is a 
noise with zero mean and its time correlation function 
is (f(t)f(t')} = fjex P (-2\t-t'\/T). We generate di- 
chotomic noise f(t) using the algorithm described in [29j | . 
In the initial equilibration time, f(t) is given with either 
+/d or —fd with an equal probability 1/2, and in a small 
time-step At, f(t) can flip to the other value with the 
probability i — i exp(— 2Ai/r), which makes the mean 
flipping time be r. 

Figure [6] A shows the normalized looping times, 
T/v(/, t)/7at(/), as a functions of r for the chain length 
iV=12 in the presence of fluctuating tension f{t) added 
to the static tension /. We consider a value of di- 
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FIG. 6: (A) The normalized looping time 7iv(/, r )/7~N(f) in 
the presence of dichotomic tension as well as static tensions 
for the chain of N = 12. The dichotomic force amplitude 
is fd = 61 fN, and the static tensions amplitudes are f — 81 
(circle) and / = 161 fN (rectangular). The looping time shows 
a minimum at optimal flipping time r r ~ 60 and ~ 300 for 
/ = 81 and 161 fN, respectively. (B) The probability P(r) 
that the dichotomic force is negative (i.e., inward direction) 
at the instant of looping as a function of r. 



chotomic noise amplitude, fa — 61 fN and two values 
of static tensions amplitudes, / = 81 (circles) and 161 fN 
(squares), which are similar to those used in a recent ex- 
periment (l2l |. For very short r, the dichotomic forces are 
averaged out, so the looping times converge to the values 
without dichotomic force. They gradually decrease with 
r until the minimum at r r 60 and « 300 for / = 81 and 
161 fN, respectively. For very long r, the dichotomic force 
rarely changes in a typical looping time, so the looping 
time goes to the average of the looping time with tension 
(/ + fa) and the looping time with tension (/ — fd). The 
average is dominated by the looping time with (/ + fd), 
so that the looping time sharply rises with r. 

Related to this, we study the probability P(t) that 
the dichotomic tension f(t) is negative (i.e., is in inward 
direction) at the instant of looping as a function of r 
(Fig. [5]£>). It has a maximum at r w 100 and r ~ 300 
for / = 81 and 161 fN, respectively, which means that, 
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FIG. 7: The normalized looping time Tn{t)/Tn in the pres- 
ence of dichotomic tension for the chain of N = 12. The di- 
chotomic force amplitudes are fd = 61 (circle) and fd = 121 
fN (rectangular). The looping time shows a minimum at op- 
timal flipping time r r ~ 40. 



for t near the t>, most of looping occurs when f(t) is 
in inward direction, and the looping time becomes the 
minimum. In the limits of very short or very long r, 
f(t) is positive or negative with equal probabilities at the 
moment of looping. The maximum of P(t) is larger for 
larger /. Evidently the minimum looping time is closely 
associated with the maximum of the P(t). 

The minimum of T is found to occur when the flipping 
time r is comparable to the diffusion time of a Brownian 
particle in the free energy F(r), r r — uL 2 / D, where a is 
a numerical value of order unity that increases with the 
tension / (30j . This phenomenon is an extension of the 
resonant activation (RA) originally found in the single 
Brownian particle crossing over a fluctuating barrier [301 ] . 
Indeed, we can regard the polymer looping in the pres- 
ence of dichotomic tension as the process of a Brownian 
particle crossing over a fluctuating free energy barrier. 

To study how the dichotomic force affects the RA phe- 
nomena, we consider different value of fd- Figure [7] shows 
the normalized looping time with fd = 61 (circles), 121 
fN (squares) in the absence of static tension (/ = 0). 
The looping time has resonant minimum at the optimum 
flipping time r r w 40 which is smaller than the case with 
static tension /. At the optimum flipping time, the loop- 
ing time 7at(0,t) is smaller for larger fa, while, at very 
long t, 7Jv(0,t) is larger for larger fa. Therefore at cer- 
tain t in between, there is a crossing point. 

In contrast to the RA of single particle, the fluctuating 
barrier heights and thus the looping time depends sensi- 
tively on the chain length N. We obtain the normalized 
looping time T/v(/, r )/T/v(/) for different chain lengths L 
in the presence of dichotomic force of amplitude fa = 61 
fN and static tension / = 81 fN (Fig. El). While the op- 
timal flipping time r r increases with L as implied by the 
relation r r = aL 2 /D, the minimum of T/v(/, T )/T/v(/) 
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FIG. 8: The normalized looping time 7jv(/, t)/Tn{J) in the 
presence of dichotomic tension as well as static tensions for 
the chains of N = 12 (circle), 18 (rectangular), and 24 (tri- 
angle). The dichotomic force amplitude is fd = 61 fN, and 
the static tension amplitude is / = 81 fN. As TV increases, the 
optimal flipping times i> increase and the minimum values of 
normalized looping time 7jv(/, T r )/7~N(f) decrease. 

decreases with L. This is because, similar to the static 
tension case, the fluctuating barrier height will be a func- 
tion of (3fdL, so that the longer chain tends to be more 
susceptible to the tension fd, and have the lower relative 
looping time. 

A recent experiment [l2| has shown that a fluctuat- 
ing tension on DNA greatly reduces the looping time for 
small r, in consistency with our results. They attributed 
the phenomenon to an increase of the effective tempera- 
ture. This may be reasonable for the small r, where the 
dichotomic noise adds to the thermal noise, but could 
not lead to the nonmonotonic resonant behavior emerg- 
ing over the entire range of r. 



IV. CONCLUSIONS 

We have studied the effects of static or time-dependent 
tension on the semiflcxiblc polymer looping using Brow- 
nian dynamics simulation. For the case of static tension, 
we have found that a minute tension as small as 100 fN 
can dramatically change the looping time, especially for 
long chains. This sensitivity is a consequence of the co- 
operativity of the chain arising from chain connectivity. 
For the case of time-dependent tension, we considered 
dichotomically fluctuating tension, where the tension in 
average changes its sign in time r. The looping time 
has a resonant minimum at an optimum flipping time 
r r , which is a nontrivial extension of the resonant ac- 
tivation (RA) of single Brownian particle. The effect of 
time-dependent tension is also more significant for longer 
chains. Our results are consistent with recent experi- 
ments for both static [ll| and a fluctuating tension [l2[ 
cases. Although we neglect the details of chain confor- 
mation outside of the loop and orientation constraints, 
etc., our model could be a basic step to understand the 
loop formation process in vivo where biopolymers are 
constantly subject to forces. This study suggests a possi- 
bility that a biopolymer can self-organize optimally uti- 
lizing the ambient nonequilibrium fluctuations. Further 
experiment using dichotomically fluctuating tension with 
various r is called for to establish the RA phenomena in 
polymer loop formation. 
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